The transverse breathing mode of an elongated Bose-Einstein condensate 
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We study experimentally the transverse monopole mode of an elongated rubidium condensate. 
Due to the scaling invariance of the non-linear Schrodinger (Gross-Pitaevski) equation, the oscilla- 
tion is monochromatic and sinusoidal at short times, even under strong excitation. For ultra-low 
temperatures, the quality factor Q = uq/to can exceed 2000, where uiq and 70 are the mode angu- 
lar frequency and damping rate. This value is much larger than any previously reported for other 
eigenmodes of a condensate. We also present the temperature variation of ujq and 70. 
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Collective excitations constitute one of the main 
sources of information for understanding the physics of 
many-body systems. Since the discovery of Bose-Einstein 
condensation (BEC) in weakly interacting bosonic gases 
10 much theoretical and experimental effort has been de- 
voted to the study of the low energy excitations of these 
systems [|[ ^, ^, |[ ||. At ultra-low temperatures, the 
measured frequencies of these excitations are found to be 
in very good agreement with the ones derived from a sim- 
ple linearization of the Gross-Pitaevskii equation, which 
describes the evolution of the condensate wave function. 

The temperature dependence of the damping and fre- 
quency of these oscillations has attracted strong atten- 



tion of both experimental 
0,|l|,|l,[l|00|ll,f 




and theoretical groups 
iJI^. Indeed, a dilute 
BEC is a very convenient system for the study of dissipa- 
tion in a macroscopic quantum object. At present, there 
still exist significant discrepancies between some avail- 
able experimental results and the corresponding theoret- 
ical predictions (see e.g. [|l5|, 0). Therefore new mea- 
surements on simple and well identified modes are clearly 
needed to test the various theoretical approaches. 

It has been known for some time that in an isotropic 
trap, the monopole (or breathing) oscillation mode may 
exhibit very unusual features. For instance, Boltzmann 
noticed that, within the framework of the equation which 
bears his name, the breathing mode of a set of particles 
confined in a 3D isotropic harmonic trap is undamped, 
irrespective of the ratio between the average collision rate 
and the trap frequency ^ . More recently it was 

pointed out that in a 2D isotropic harmonic trap, the 
breathing mode of identical bosons interacting via a con- 
tact potential is also undamped This property, 
which holds at any temperature, can be related to the 
existence of a hidden symmetry of the problem described 
by the two-dimension Lorentz group S0(2,l) |25| . 

The purpose of the present work is to investigate ex- 
perimentally the transverse breathing mode of an elon- 
gated condensate. We find that this mode exhibits very 
unique features, such as an extremely large quality factor 
Q — luq /70 ^ 2000 for temperatures much lower than the 
condensation temperature. Here uq and 70 are the angu- 



lar frequency and damping rate, respectively. This qual- 
ity factor, much larger than any previously reported value 
for an oscillation mode of a gaseous BEC, should provide 
a sensitive test of the various theoretical models of dissi- 
pation at ultra-low temperature. Another very unusual 
property lies in the fact that the measured frequency loq 
is nearly independent of temperature. Furthermore, we 
find that the transverse area of the condensate undergoes 
a pure, isochronal oscillation even for strong excitation. 
This striking behavior emerges from a scaling invariance 
of the time-dependent Gross-Pitaevski equation. 

Our ^^Rb condensate is formed by radio- frequency (rf) 
evaporation of 10^ atoms in a loffe-Pritchard magnetic 
trap. The initial temperature of the cloud pre-cooled us- 
ing optical molasses is 100 fiK. We vary the final tempera- 
ture of the gas and the number of atoms in the condensate 
by adjusting the final rf, j/f . The condensation threshold 
is reached at = 290 nK, with = 8.4 10^ atoms [|6|. 
The coldest sample reliably reproducible {T ^ 40 nK) is 
a quasi-pure condensate with Nq = 10^ atoms. It is ob- 
tained using h'f = vq + ^ kHz, where vq is the rf at which 
the trap is emptied. 

The trap has a longitudinal frequency uJz/(2tt) — 
11.8 (1) Hz for the atoms prepared in the stretched state 
TO = 2 of the 5 S'1/2, F = 2 ground state. The transverse 
magnetic gradient is b'^ =^ 1.33 Tm^^ and the transverse 
frequency u!± is adjusted by varying the bias field Bq at 
the center of the trap. We choose Bq = 0.844 10"'' T, 
which leads to u;j_/(27r) = 182.6 (4) Hz Both fre- 

quencies ujz and uj± are measured by monitoring the os- 
cillation of the center-of-mass of the atom cloud. 

Once the condensate is formed, we excite the trans- 
verse breathing mode by modifying the bias field from 
-Bo to Bq, which changes the frequency uj± to uj'j_ > ll!±. 
After a time r ^ 2tt /uj'^, we set the bias field back to 
its original value. We let the cloud oscillate freely in the 
magnetic trap for an adjustable time t and then measure 
the transverse density profile of the condensate after a 
period of free expansion. In this pursuit, we suddenly 
switch off the magnetic field, allow for a 25 ms free-fall, 
and image the absorption of a resonant laser by the ex- 
panded cloud. The imaging beam propagates along the z 
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FIG. 1: Variation of A{t) = Rl{t) + Rl{t), after a strong 
excitation of the transverse breathing mode. The line is a 
sinusoidal fit of the data (ijj'^I2-k = 250 Hz, r = 800 /is). 

axis. We fit the density profile of the condensate assum- 
ing the parabolic shape of the Thomas-Fermi distribution 
and extract the radii and Ry in the plane z = 0. The 
temperature is obtained from a Gaussian fit of the un- 
condensed part of the distribution. 

Within the Thomas-Fermi approximation, the evolu- 
tion of the trapped condensate after the excitation of the 
transverse breathing mode is well described by a time- 
dependent scaling transformation |2^, The spa- 
tial density p(r, t) at time t is deduced from the spatial 
density at time (i.e. just after the transverse frequency 
is set back to its value LO\}j by: 
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The evolution of the scaling parameters A_l and A^ is 
given by H H: 



A_L — —3-^ I A I 



A I Aj. 



We start our analysis with the evolution for relatively 
short times t <C lixjiOz- In this case one can neglect the 
variation of \z and integrate the equation for Aj^: 

1 /2 

A_L(t) = (cosh(a) + sinh(a) sin(2a;^< — /3)) ' , 

where the parameters a and /3 depend on the initial con- 
ditions. One remarkable consequence of this solution is 
that the quantity A{t) = Rl{t) + Rl{t) = \\{t) A{0) al- 
ways undergoes a sinusoidal oscillation at frequency 2lj^^ 
irrespective of the strength of the excitation . This 
is confirmed by the results of fig. ^ where we observe 
a sinusoidal oscillation of A{t) with a factor 3 between 
the maxima and the minima. The frequency of the os- 
cillation is indeed found to be equal to 2wj^, within our 
experimental uncertainty (±0.5 %). 

At longer times one has, in principle, to take into ac- 
count the variation of Xz- We restrict our analysis to 
small amplitude oscillations and write Xi = 1 -f , where 
£±56^ 1. The linearized equations of motion 

+ ojI (4ex + e.) = h + ujI {2e^ -f 3e,) = 



give rise to two modes H, ^. One is a fast mode of 
frequency 2uj± corresponding to e±/ez = ujj_/uJz. The 
other one is a slow mode of frequency y^5/2 lOz with 
t^jtz ~ —1/4. Our excitation scheme corresponds to 
ej^ 7^ at time while the three other quantities e^, 
tz and tz are approximately zero at initial time. Under 
these conditions, we mainly excite the fast mode 2uj^ 
since the relative weight of the slow mode in the evolution 
of A_L is only Wz/(\/40tjj_) < 0.01. For this fast mode, the 
condensate dynamics consists essentially in a transverse 
monopole oscillation since Ez ^ e^. 

We investigate the oscillation and the damping of this 
mode on a long time scale (up to 0.7 s). An exam- 
ple is given in fig. | for T = 40 (±20) nK (note 
that for such cold clouds, corresponding to a quasi-pure 
condensate, the temperature can only be inferred from 
the final rf used for the evaporation). The evolution 
of A(t) is well fitted by a damped sinusoidal function, 
Ao + Mo cos(woi + 0) e'T«*, with Mo Mo = 0.063 (4), 
tJo/27r = 366.3 (5) Hz and 70 = 1.2 (2) s^^. This corre- 
sponds to a quality factor Q = wq /70 ^ 2000 much larger 
than any previously reported for other eigenmodes of a 
BEG. For instance, measurements performed with a TOP 
trap [lOz = 2\/2 wj,) gave Q ~ 200 for the lowest to = 
mode 0. For a cigar shape trap, measurements of the 
low frequency mode discussed above (w = •\/5/2 lOz) led 
to Q ~ 80 [H. A smaller value (Q ^ 25) was measured 
for the scissors mode |9|. The present quality factor is 
also one order of magnitude larger than what we find for 
the transverse quadrupole mode under the same experi- 
mental conditions ||3|] . 

We measure the oscillation frequency and the damp- 
ing rate over a large temperature range compatible with 
the detection of a condensate. The results are given in 
fig. ^, together with the total number of atoms N . The 
measurement of 70 is made from the decay of the oscil- 
lation amplitude between and 200 ms. For the points 
of lowest and highest temperature in figs. 3b and 3c, we 
calculate, including mean field and finite size corrections 
T/Tc = 0.3 and 0.8, respectively. 

The frequency is found to be very close to the 
Thomas- Fermi prediction 2ci;j_, and it varies only slightly 
with temperature (fig. 3b). In particular, we do not ob- 



serve the 10 % relative decrease predicted in |14 when 
T/Tc increases from to 0.8. 

The damping rate 70 varies approximately linearly 
with temperature in the range 50 — 200 nK. The linear 
variation of 70 for k^T larger than hw^ (^ 20 nK) and 
than the chemical potential (~ 60 nK for the coldest 
point in fig. ^ is expected if one assumes either Beliaev- 
type or Landau-type damping of the monopole excitation 
(see e.g. [^). In the Beliaev mechanism, the monopole 
excitation decays into two excitations with a lower en- 
ergy. For Landau damping, it annihilates in a collision 
with a thermal excitation to give rise to another exci- 
tation. The calculation of the corresponding damping 
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FIG. 2: Long time evolution of A(i)IA at low temperature 
[0,9 ms], [200 ms,209 ms] ... The oscillation is well fitted by a 



(40 nK), where A is the average of A{t) over each time interval 
damped sinusoid. Here uj'j_/2ti = 230 Hz and r = 75 /xs. 
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FIG. 3; (a) Total number of atoms A'^, (b) breathing mode 
frequency ujo/uj± and (c) damping rate 70, as a function of 
temperature. Each point in (b) and (c) corresponds to i/t — 
vo=Q kHz, 9 kHz, . . . , 33 kHz = 230 Hz, r = 150 

rate is quite complex for inhomogeneous systems and can 
lead to drastically different results for modes with similar 
energies, depending on the stochastic or regular nature 
of the relevant thermal excitations . In this respect 
it is worth noting that the frequency of this transverse 
breathing mode is an integer multiple of the trap oscil- 
lation frequency, which is usually not the case for other 
modes of the condensate. 

For a very low temperature (T ~ 40 nK) and an in- 
creasing excitation strength, we observe a significant de- 
viation of the long term behavior of A{t) with respect 
to the simple damped sinusoid shown above. Fig. ^ 
shows the oscillation amplitude for the weak excitation 
(•, T = 75 /is) together with two sets of data obtained for 
Ti = 300 /iS (□) and T2 — 450 /is (a). One clearly sees 
that for strong excitation the amplitude of the oscilla- 
tions passes through a local minimum around t — 300 ms 
before increasing again, and decaying finally after 700 ms. 
We fit our experimental data with the phenomenological 
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FIG. 4: Long term evolution of the amplitude of the breath- 
ing mode for different excitation strengths {•■.t = 75 /is, □: 
Ti = 300 /xs, k:T2 — 450 /xs). Each point ai t — 0, 100 ms, 
. . . , 700 ms represents the amplitude of the breathing mode 
measured during the period [t,t + 8 ms]. The fit function is 
indicated in the text. Here lj'^/2tv — 230 Hz. 

function A{t) ^ Aq + 5A{t) cos{ujQt + (/>) where the am- 
plitude dA{t) is now slowly modulated in time, in addi- 
tion to the exponential decay already discussed: 6A{t) = 
SAq (1 — sin^ {5 Lu t)) e~^*. The best fit parameters are 
for these two cases 71 = 1.48 (15) s"-^, (5wi/(27r) ~ 
1.08 (4) Hz, ^1 = 0.5(1), and 72 = 1.35 (15) s^^, 
(5w2/(27r) = 0.96 (4) Hz, 6 = 0.6 (1). 

Two scenarios can be invoked to explain this evolution 
of A{t), which is not predicted by (|l|). A first expla- 
nation could consist in assuming that we actually ex- 
cite two modes whose frequencies ~ loq differ by a small 
amount and that we observe the beating between these 
modes. However, if we assume that these two modes 
are excited linearly, the relative amplitude of the beating 
should not depend on the excitation strength, contrary to 
the observed. Furthermore, this explanation is not con- 
firmed by a linear analysis of the time dependent Gross- 
Pitaevskii equation for our experimental conditions [^ . 
A second, more plausible, scenario involves a non-linear 
process, in which the monopole excitation may be con- 
verted in a reversible way into a sum of other excitations 
of frequency loj. This process is efficient if it is resonant: 

^0 ^ E, 



In this context, the modulation in fig. 



can be interpreted as a parametric oscillation. 

As a possible candidate for such a parametric (or Be- 
liaev type) oscillation |3^, one may consider the excita- 
tion of the pair of modes I = 2,m = ±1 P,i, p|. The 
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frequency of each mode is (w^ + tj^Y^^ so that the reso- 
nance condition given above is indeed nearly fulfilled. By 
imaging the condensate along the y direction, we observe 
a weak excitation of this mode. However, no increase 
in its amplitude is found around t =300 ms when the 
amplitude of the breathing mode collapses (cf. fig. ||). 

Alternatively one may consider the excitations con- 
sisting in a pair of phonon-like excitations propagating 
along z with momentum k and — fc, each with an en- 
ergy HuJo/2 ~ Hcl!± To test this latter hypothe- 
sis, we perform the same experiment for various values 
of uj±, so that it spans the whole interval between two 
successive phonon modes (typically 6 Hz in this energy 
domain). However, expecting a resonance when wq/S co- 
incides with the frequency of a phonon mode, we do not 
observe any significant change in the evolution of A{t). 

In any case we note that for the long oscillation times 
involved here, small non-linear terms in the confining po- 
tential (varying as x'^z and y'^z) may play a significant 
role. Further work investigating this possibility is under 
progress in our laboratory. 

In summary we have presented in this letter a study 
of the transverse breathing mode of an elongated Bose- 
Einstein condensate. For low temperatures and short 
evolution time, the observed isochronal and sinusoidal 
oscillation is in excellent agreement with the prediction 
based upon the scaling invariance of the Gross-Pitaevski 
equation. However, the temperature dependence of the 
frequency and damping of the oscillation reveals, on 
the contrary, some unexpected features. First, the ob- 
served damping is comparatively ten times smaller than 
what was previously reported for other modes of a BEG 
(quality factor of 2000 instead of 100-200). Second, the 
frequency of the breathing mode is found to be quasi- 
independent of temperature, in contradiction with recent 
theoretical predictions. Finally the long time evolution 
of this mode when strongly excited reveals an unexpected 
non-linear mixing phenomenon. This could be explored 
further using the formalism of |^ , by looking at a possi- 
ble dynamical instability of the Gross-Pitaevskii equation 
in this regime of strong excitation. 
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